Role of fluctuations in membrane models: thermal versus non-thermal 
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We study the comparative importance of thermal to non-thermal fluctuations for membrane-based 
models in the linear regime. Our results, both in 1+1 and 2+1 dimensions, suggest that non- 
thermal fluctuations dominate thermal ones only when the relaxation time r is large. For moderate 
to small values of r, the dynamics is defined by a competition between these two forces. The 
results are expected to act as a quantitative benchmark for biological modeling in systems involving 
cytoskeletal and other non-thermal fluctuations. 
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INTRODUCTION 

A long standing tradition in the field of non- 
equilibrium statistical mechanics, especially those 
in connection with biophysical dynamics [IH1] has 
been to look at the high temperature limit of all 
such models where thermal fluctuations dominate 
every other stochastic mode in the system. Such an 
approach generally has the advantage of having to 
deal only with a thermal noise that dominates all 
non-thermal motions in the mesoscopic limit. Both 
theoretically [THI] and experimentally [5--8J, the 
subject is replete with such examples, and often 
with agreeable results too, where such an assump- 
tion agrees with the physiology of the real system. 
Problems, though, surface when the approach con- 
tradicts its basic tenet in that the system has a 
non-trivial non-thermal force that competes, and 
sometimes even out-competes, the thermal fluctu- 
ations at the experimental setting [5HT2]. In the 
present article, our focus is on the latter situation 
when the system has a non-trivial component of 
non-thermal fluctuation eg. cytoskeletal fluctua- 
tion [SHTT] whose contribution can not be neglected 
in comparison to its thermal counterpart. 

A similar, if not identical, question in connection 
is whether or not there exists a non-equilibrium 
version of the fluctuation-dissipation (FD) theo- 
rem in the mould of an Jarzynski relation . 
Our approach, though, will rely on a relative 
comparison of strengths between thermal and 
non-thermal (stochastic) fluctuations connected 
to identical reservoirs. 

A remarkable range of such non-thermal fluctu- 
ations are visible in the biological world, a typ- 
ical case in hand being the fluctuations due to 
bio-chemical and molecular movements, popularly 
referred to as 'cytoskeletal fluctuations'. Recent 
studies [16] have traced the origin of such (gener- 
ally) non-thermal cytoskeletal dynamics to a com- 



bination of cell crawling, adhesion between lipid 
bilayers and randomness in the responses of molec- 
ular motor proteins. In this article, though, we 
would not pay so much of attention as to the origin 
of such motions; rather our focus would be on a rel- 
ative comparison of strengths between a Brownian 
thermal noise against a randomness that is mostly 
non-thermal in origin eg. cytoskeletal fluctua- 
tions and whose strength is proportional to y/T e g, 
T e g being a 'non-equilibrium temperature'. The 
Langevin and Fokker-Planck type models based on 
such a description [2J [31 [T7] generally utilize the 
fact that the amplitude of thermal motion thor- 
oughly dominates all other mesoscopic level fluctu- 
ations that are otherwise present in the dynamics. 
Although such an approach works remarkably well 
in many cases [5} [B] , and we are mostly focusing on 
native biological systems involving cytoskeletons, 
serious complications arise in cases where the cy- 
toskeletal fluctuations are non-negligible. The sim- 
ple question one would like to ask here is the fol- 
lowing: what are those situations when such a non- 
thermal motion competes with a thermal fluctua- 
tion attached to a heat bath 18J? A notable contri- 
bution to the understanding of such non-thermally 
initiated fluctuations owes to Prost and Bruinsma 
|19j . They showed that in systems involving ac- 
tive membranes the long-wavelength limit of the 
fluctuation spectrum follows a power law behavior 
that is quite different from a membrane in ther- 
mal equilibrium. Using this approach they were 
successful in analyzing the experimental results on 
red-blood cells [T2] where cytoskeletal fluctuations 
were convincingly seen to be dominating the ther- 
mal ones. Similar studies dealing with the effect of 
non-thermal motions in the intra-cellular traffick- 
ing involving membrane fluctuations [20j , or in the 
travelling waves observed due to protein activity in 
a flexible membrane [ST], as well as in micropipet 
experiments using an activated membrane surface 
lend credence to the belief that active non- 
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thermal motion can indeed have a dramatic effect 
on the overall dynamics of a non-equilibrium sys- 
tem, especially when there is a competition with a 
fluctuations arising from a stochastic thermal bath. 

In this article, we study the regime when the 
two noise modes, thermal and non-thermal, com- 
pete with each other. Ramifications of such 
comptetitions are amply evident in experiments 
[9HI21 [22) [23] and our attention here is to quan- 
tify the limits of the parameter that defines such 
a regime. As we would later see, the parameter in 
question is the relaxation time r that, in a way, 
offers a complementary description to that studied 
in (T5]. We study the model in both asymptotic 
limits of t and predict values (indeed limits) of r 
for which thermal noise dominates (or competes 
with) its non-thermal counterpart. 
Our basic model consists of a large flat d+1 
dimensional membrane driven by thermal fluctu- 
ations where the free energy is Helfrich-like |24) . 
involving a term that is proportional to the energy 
cost due to surface tension and another term that 
represents the bending energy associated with 
the existence of curvature in a membrane. In 

our nomenclature, we will represent the surface 
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tension term as F sur f = ^ J d d r(VZ) and a 
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curvature term Fb en d = / d d r(\/ 2 Z) . However, 
in the present article, we focus only on the 
linear regime of such an intrinsically non-linear 
potential (in the mould of [17). As can be easily 
seen, this linear regime can be arrived at as 
a high-temperature expansion of the Kramer's 
type two-state potential studied in [TJ |2]. As 
is always true of linear models, ours too would 
be constrained by the fact that exact measure- 
ments of amplitudes of perturbations as well 
as enacting proper boundary conditions would 
suffer. However, the gross aim being a qualitative 
understanding of the role of noise (both thermal 
and non-thermal) in perturbed membranes, we 
forego this without any effective loss of generality. 

Our interest is in the activated dynamics of this 
membrane in the presence of a thermal, together 
with a non-thermally fluctuating noise. We de- 
fine the thermal part of the noise using the stan- 
dard FD theorem, while the non-thermal part in- 
cludes a time-decaying component. The finite de- 
cay in the non-thermal noise distribution simply 
indicates the presence of a finite relaxation time 
r as opposed to a an infinite relaxation time for 
a white noise source. One might consider more 
complicated noise sources as well, like a non-local 
spatial correlation in its distribution. We, how- 



ever, stick to the minimalist model since a spatial 
non-locality does not change the qualitative out- 
come, apart from complicating the algebra that is. 
In the following description, the thermal noise is 
represented by r)^(x,t) and the non-thermal noise 
by 77nth(£, t)- The over-damped dynamics of such a 
membrane (we consider both 1+1 and 2+1 dimen- 
sional systems individually to illustrate the gen- 
erality of the argument) defined by a stiffness B, 
surface tension T and mobility 7 is given by 



1&Z %^ = -£V 4 z(x,;) + TV 2 z(f,<) (i) 

+ nonlinear terms + -F oxt 

.Fext is the space-time dependent external force im- 
pressed on this system that, in principle, can be 
stochastic. Defined generally, we define the force 
as a monotonically decaying function of time as 
follows 



< F ext (x,t)F ext (x f ,t f ) >= F cxp[- 



\l—Q]6(2-3?) 

T 

(2) 

Our objective is to study the spatio-temporal prop- 
erties of this model both in 1+1 and 2+1 di- 
mensions for the two asymptotic regimes r — > 
and r — > 00 and compare such responses with a 
stochastically (thermal) driven string/membrane. 



NON-THERMAL FLUCTUATIONS IN A 
1+1 DIMENSIONAL STRING 

We start from a linearized version of a (1+1 dimen- 
sional) membrane model driven by a combination 
of an external time dependent non-thermal noise 
77 n th and a Gaussian thermal noise 77th ■ 



l dZ{ Q^ = -BV 4 Z($,t)+TV 2 Z{$,t)+r) nth +r) th 

(3) 

where the noises 77 n th and 77th are defined as follows: 



< Vnth(x,t)r/ nth (x ,t') > 



F 8(x* -xty) 

D 6{t - t') 5{x - x') 




Since the thermal and non-thermal noises are 
uncoupled [35] (although not mutually exclusive; 
in fact the main aim of our study is to ascertain 
the dominant mode when both are simultaneously 
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present), we will decouple equation (|3| in to sepa- 
rate thermal and non-thermal modes without any 
loss of generality. Note that this is possible only 
when the non-thermal amplitude is completely 
independent of ('non-equilibrium') temperature 
as is the case in most biological systems. For 
exact quantisation of the numbers to be evalu- 
ated henceforth, we would use parameter values 
as in pi]: B=11.8 k B T, T=5650 k B T /im~ 2 , 
7 = 4.7 x 10 6 k B T s fim-' 1 and D = 2k B T M. 

Our working model for the non-thermal case would 
then be 



1 dZ^A = _ B ^ Z{ ^ t) + tV 2 Z(x, t) + Vnth 
at 

(5) 

where F ~ 1/r due to normalisation constraint. 
Solving the above equation in the wave vector 
space, we evaluate the two-point temporal corre- 
lation function as 



< ZAt)Z M) >= § 

T[a(k) - ±] 



(6) 



{ 



cxp[-(^)] 1 
a(k) + i 2a{k) 



exp[-a(fc)(i' - <)]} 



for t! > t, where a(k) 



_ Bk^+Tk 2 
7 



This gives us 



the detailed two-point structure function in 1+1 
dimensions as 



7 -5/4 T -l/4_gl/4 



4T 



< Z nth (x,t)Z nth (x,t') >- 
1 1 , 
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-5/4 T -l/4 S l/4 



X ( 



i exp[ 



4T 
(f-t). 



(8) 



i r°° dk e -«(£)i*'-*i 



where in the above equation we have replaced Z 
by Z nt h to indicate that it is the auto-correlation 
function for the non-thermal scenario. 

We now move on to a similar looking model as 
defined in cqn (4) but with a thermal noise rjth 
now replacing the non-thermal noise n n th' 



7 ^^^ = -BV 4 Z ti (£,t)+TV 2 Z th (f,t)+r, th (x,t) 

(9) 

where 77th is as defined in eqn (pj). Proceeding 
as before, the two point structure function for 
the thermally perturbed membrane in the Fourier 
transformed k-space reads as 



< Z th (k,t)Z th (-k,t') >-- 



D Q e -«WI*'-*l 



r 



2a(k) 



(10) 



In the limit X — > and \t' — t\ —> 00, that is for 
very close spatial points evolved until the system 
reaches the stationarity limit, we get 



< Z th (x,t)Z t h(x,t') >= 



D f°° dkexp[-a(k)(f -t)] 



a(k) 
(11) 

^FCquation (11) above resembles the second part 



(t' - t) exp[-Xy/-/j, + \ exp[~Xy^JL 
x exp[ — \{- 



of equation (8) (eqn (11) approaches eqn (8) 



1 f 00 dk e ikX - a (k)\t'-t\ 



V Jko 2^ a (k) 



] in the limit r — > 0) which implies that it is 
'essentially the first part of the same equation 
(8) that will decide which noise dominates the 
system when both are simultaneously present. 
We would revisit this question soon after solving 
cqn ( 11 ) in the stationary state limit (\t'— t\ — > 00). 



where /x ± = |[-^ ± \Jj^ -4], with e = j/t 

and X = {-g^) 1 ^ 4 ^ — x'\. kg relates to the smallest 
length scale in the system, which, in effect, turns 
out to be the lattice size A (fc = 

The above equation specifies the two-point tempo- 
ral auto-correlation function defined at the same 
spatial point (X=0) as 



_ foo dk cxp[q(fc)(t'-t)] 
JO 2tt 

be represented as the following integral equation 



Let us define J = f n °° This can 

■'" a(k) 



dJ 



= [ *p-«(fe)|t'-«l 



(12) 



d\t'-t\ J 2vr' 

which can be approximately solved in the limit 
fco — > to get 
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( ns) relaxation times are not unknown [53] . In the 
limit of smaller relaxation times, that is transient 
perturbations, the non-thermal (eg. cytoskeletal) 
fluctuations compete against all thermal modes. 
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NON-THERMAL FLUCTUATIONS IN A 
2+1 DIMENSIONAL MEMBRANE 



Taking cues from the linearized model for non- 
thermal fluctuations eqn ^ and its eventual k- 
space equivalent, we arrive at 



FIG. 1: The outset shows the variation of J(\t' — t\) 
with varying relaxation time r for \t' — 1\ ~ r while the 
inset shows the ratio of the non-thermal to thermal 
contribution R(t). 



< Zt(t)Z_M') >-- 



Fa 

7 2 Kfc)-i] 



(15) 



x { GXP L ( X )] - — cxp[-a(fc)(f' - t)}} 
a(fc) + i 2a(k) 



dJ(\t'-t\) 
d\t'-t\ 



T 2 



where K(n, x) represents the n th order Bessel func- 
tion of the second kind. In the \t' — 1\ — > large limit, 
the above expression reduces to 



J(\t'-t\) 



Erf(- 



T 



B\t'-t\ '78^7 



(14) 

Erf(x) represents the incomplete Gamma function 
(also known as the error function) that is repre- 
sented as Erf(x) = ^= f Q du e~ u . 
The comparative strength of the non-thermal 
versus the thermal noise is encapsulated in 
the ratio R = <Z "% ( '^f; th( '' O> as could 

be evaluated from equations (8) and (13). 
With the parameter values previously specified, 
^J{\t' -t\)\ w _ t \^ T ~ 10- 12 , implying R > 
(inset of the figure). As can be clearly seen from 
Fig. [T] (reminiscent of experimental observations 
too [23 ), in the limit of large r, J becomes 
independent of the relaxation time while R too 
starts to saturate. At this point (r — > large), 
non-thermal noise starts taking over the thermal 
fluctuations. With the parameter values used, r 
is of the order of a few seconds (as shown in the 
figure), although much larger ( hours) or smaller 



As in the previous section, we utilize this infor- 
mation to evaluate the 2+1 dimensional structure 
function 



< Z(x,t)Z(^,t') > 
d 2 k 



(16) 



AkX cos 



(2nY 



1 



47TT7 2 

1 

2a(k) 



dkkJ {kX){ 



a(k) + i 



-a(k)(t'-t) 



} 



where Jo(kX) is the zeroth order Bessel function 
of the first kind. Resorting to the usual limit 
A^0, we get J (fcA) ^ 1 - ^ +0(fc 4 A 4 ). As 
can be easily seen, eqn p7] ) is perfectly integrable 
for the first term in this equation while the second 
term shows a mild logarithmic divergence that 
can be taken care of by assuming a minimum 
cut-off length scale for the system, a length that 
simply represents the lattice size of the discretized 
system. The basic conclusion remains unchanged 
though, that for large values of the relaxation 
time r the non-thermal fluctuations will dominate 
the thermal ones while for small values of r 
the dynamics will be decided by a competition 
between these two. 
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CONCLUSIONS 

In this paper, we have discussed the relative 
importance of thermal noise with respect to 
non-thermal fluctuations. The results give us a 
clear quantitative basis for considering or neglect- 
ing non-thermal fluctuations in representative 
biological models. This is of utmost importance, 
since a traditional disturbing tendency in such 
studies has often been to neglect non-thermal 
modes in favor of thermal noise. Although, this 
entails the disconcerting possibility of having 
to consider contributions from an often less 
understood source of fluctuation, like cytoskeletal 
fluctuations in the analyzes of microtubule related 
studies, an advantage on the hindsight is that of a 
theoretical clarity. A knowledge of the relaxation 
mechanism that hopefully would be available from 
the experimental studies would be a huge boost to 
the complementary theoretical modeling of such 
systems. Work is presently underway to evaluate 
quantitative modifications in the presence of 
competing non-linearities in the dynamics. 
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